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A LOCKING-FREE WEAK GALERKIN FINITE ELEMENT METHOD 
FOR ELASTICITY PROBLEMS IN THE PRIMAL FORMULATION 

CHUNMEI WANG*, JUNPING WANGt, RUISHU WANG*, AND RAN ZHANG§ 


Abstract. This paper presents an arbitrary order locking-free numerical scheme for linear 
elasticity on general polygonal/polyhedral partitions by using weak Galerkin (WG) finite element 
methods. Like other WG methods, the key idea for the linear elasticity is to introduce discrete weak 
strain and stress tensors which are defined and computed by solving inexpensive local problems 
on each element. Such local problems are derived from weak formulations of the corresponding 
differential operators through integration by parts. Locking-free error estimates of optimal order are 
derived in a discrete iL^-norm and the usual L^-norm for the approximate displacement when the 
exact solution is smooth. Numerical results are presented to demonstrate the efficiency, accuracy, 
and the locking-free property of the weak Galerkin finite element method. 
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1. Introduction. In this paper, we are concerned with the development of effi¬ 
cient new numerical methods for linear elasticity equations by using the weak Galerkin 
finite element method recently developed in [lIJIllllllIl]. Let n c (d = 2,3) be 
an open bounded and connected domain with Lipschitz continuous boundary T = dQ 
of an elastic body subject to an exterior force f and a given displacement boundary 
condition. The kinematic model of linear elasticity seeks a displacement vector field 
u satisfying 

(1.1) — V ■ tT(u) = f, in O, 

(1.2) u = u, on r, 

where (7(u) is the symmetric Cauchy stress tensor. For linear, homogeneous, and 
isotropic materials, the Cauchy stress tensor is given by 

ct(u) = 2/x£(u) + A(V • u)I, 

where e(u) = i(Vu-f-Vu^) is the linear strain tensor, /i and A are the Lame constants. 
For linear plane strain, the Lame constants are given by 

Ely E 

^ + ^^2(1+7)’ 
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where E is the elasticity modulus and v is Poisson’s ratio. 

Denote by (•, •) the L^-inner product in either L^(0), or as 

appropriate. A weak formulation for (HI) in the primal form reads as follows: Find 
u G satisfying u = u on F and 

(1.3) 2(/re(u),e(v)) + (AV ■ u, V • v) = (f,v), Vv G 

where is the usual Sobolev space defined by 

H^{n) = {u : v£ L^{VL), Vv G 

and 1(^0 (fl) is the closed subspace of consisting of all the functions with van¬ 

ishing boundary value. 

The main objective of this paper is to study an application of the weak Galerkin 
finite element method [m EH EH Eol [la [IS] to the linear elasticity problem (HU- 
(EEl) based on the primal formulation dia. The weak Galerkin (WG) refers to a 
generic finite element technique for partial differential equations where differential 
operators are approximated or reconstructed by solving inexpensive local problems 
on each element. Such local problems are often derived from weak formulations of 
the corresponding differential operators through integration by parts. Recent work on 
WG has revealed that the concept of discrete weak derivatives offers a new paradigm 
in the discretiztion of partial differential equations. The resulting numerical schemes 
often possess a great robustness in stability and convergence for which other com¬ 
peting methods are hard to achieve. For the linear elasticity problem (11.311 . we shall 
demonstrate that the WG numerical approximations are not only accurate and robust 
with respect to the polygonal/polyhedral partition of the domain, but also naturally 
“locking-free” in terms of the Lame constant A. This is a result that the standard 
conforming finite element method does not have. 

“Locking” refers to a phenomenon of numerical approximations for a certain prob¬ 
lems whose mathematical formulations involve a parameter dependency. For the linear 
elasticity problem, the parameter is the Poisson ratio ly. For i/ close to A (i.e., when 
the material is nearly incompressible), it is well known that various finite element 
schemes, such as the continuous piecewise linear elements, results in poor observed 
convergence rates in the displacements. In 1983, Vogelious [18] showed absence of 
locking for the p-version of the finite element method on smooth domains. Later on, 
Scott and Vegelious m proved that no locking results when polynomials of degree 
fc > 4 are used on triangular meshes. However, Babuska and Suri |2] found, for 
conforming methods, locking cannot be avoided on quadrilateral meshes for any poly¬ 
nomial of degree fc > 1. In the discontinuous Galerkin context, Hansbo and Larson 
m proved that the numerical approximation arising from a discontinuous Galerkin 
method is locking-free for any values of fc > 1, also see E3] for the case of fc = 1. 
In |S] , Daniele and Nicaise designed a locking-free discontinuous Galerkin method for 
composite materials featuring quasi-incompressible and compressible sections. 

Locking occurs because for the limit case ^ (or A = oo), the exact solution 
of the displacement must satisfy the constraint V • u = 0. But locking is not a 
difficult issue to resolve for the linear elasticity equation. One possibility is the use of 
mixed methods by reformulating the linear elasticity equation as a generalized Stokes 
equation (see umii] and the references cited therein). Specifically, by introducing a 
pressure variable p = AV ■ u, the elasticity problem (I1.1D - (I1.2D can be reformulated as 
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follows: Find u € [i/^(r2)]^ and p G satisfying u = u on F, the compatibility 

condition X~^pdx = Jp u ■ nds, and the following equations 

(1.4) 2(/ie(u),e(v)) + (V ■ v,p) = (f,v), Vv e 

(1.5) (V ■ u,g) - (A”V,g) = 0, VgeLg(fl). 

Here n is the outward normal direction on F, and Fo(^) i® closed subspace of L^(n) 
consisting of all functions with mean-value zero. Consequently, any finite elements 
that are stable for the Stokes problem would provide a locking-free approximation 
for the linear elasticity problem but at the cost of solving a saddle-point 

problem with an additional pressure variable. 

The weak Galerkin method can also be applied to the linear elasticity problem 
based on the mixed formulation (ll.4D - (|1.5D . In principle, such applications should 
yield locking-free numerical approximations for the displacement variable. Surpris¬ 
ingly, we found that the weak Galerkin finite element method based on the primal 
formulation (11.31) is equivalent to the weak Galerkin when applied to the mixed formu¬ 
lation (fni-lfTSl). This equivalence indicates that the weak Galerkin approximation 
arising from the primal formulation (11.31) might be locking-free. The main goal of this 
paper is to provide a rigorous mathematical justification to the locking-free nature of 
weak Galerkin. Specifically, we shall perform the following tasks in this study: (1) 
propose two weak Galerkin finite element schemes, one based on the primal formu¬ 
lation m and the other based on the mixed formulation (ll.4l) - (ll.5l) : (2) show that 
the two weak Galerkin finite element schemes are equivalent; (3) establish a locking- 
free convergence theory for the WG scheme based on (fT:il)-(IT31): and (4) numerically 
demonstrate the accuracy and the locking-free property of the proposed weak Galerkin 
finite element methods. The main mathematical challenge of this research lies in the 
error estimate for shape-regular finite element partitions consisting of arbitrary poly¬ 
gon or polyhedra. The corresponding technicality is represented by the inequality 
(19.161) in Lemma [mni which is a result of a creative use of the Korn’s inequality and 
the domain inverse inequality [21] . 

Throughout the paper, we will follow the usual notation for Sobolev spaces and 
norms [5] . For any open bounded domain D C with Lipschitz continuous boundary, 
we use II • ||s_D and | • |s,d to denote the norm and seminorms in the Sobolev space 
H^{D) for any s > 0, respectively. The inner product in H^{D) is denoted by (•, ■)s,d- 
The space H^{D) coincides with L‘^{D), for which the norm and the inner product are 
denoted by || • Ho and (•i ■)d, respectively. When D = ft, we shall drop the subscript 
D in the norm and inner product notation. 

The paper is organized as follows. Section 2 is devoted to a discussion of weak 
divergence and weak gradient for vector-valued functions as well as their discrete 
analogues. In Section 3, we present a weak Galerkin finite element method for the 
linear elasticity problem based on the primal formulation (11.31) . In Section 4, we 
describe another weak Galerkin finite element method based on the mixed formulation 
dLiii-drsi). It is also shown in Section 4 that the two weak Galerkin methods are 
equivalent. Section 5 is devoted to a discussion of stability conditions (i.e., the inf- 
sup condition and coercivity estimates) for the mixed weak Galerkin finite element 
method. In Section 6, we prepare ourselves for error estimates by deriving an identity. 
Section 7 is devoted to the establishment of an optimal order error estimate in a 
discrete iJ^-norm. In Section 8, we use the usual duality argument to derive an 
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optimal order error estimate in the L^-norm for the displacement variable. In Section 
9, we derive some supporting tools and inequalities useful for error analysis. Finally, 
in Section 10, we report some numerical results that demonstrate the accuracy and 
locking-free nature of the proposed weak Galerkin finite element methods for the linear 
elasticity problem. 

2. Weak Divergence and Weak Gradient Operators. The divergence and 
gradient operators are two primary differential operators in the variational problem 
(HU-(ESI). The goal of this section is to review the definition and computation of 
weak gradient and weak divergence operators which have been introduced and studied 
in applications to other partial differential equations [mini [22]. 

Let K he a, polygon in 2D or polyhedra in 3D with boundary dK . Define the 
space of weak vector-valued functions in K as follows 

V{K) = {v = {vo,Vb} : Vo € S [L^dK)^}, 

where Vp and Vf, represent the values of v in iF and on the boundary dK, respectively. 
Note that v;, is not necessarily related to the trace of Vq should it be well-defined. 
Denote by {■,-)dK the standard inner-product in either L'^{dK) or [L‘^{dK)]‘^ , as 
appropriate. 

Definition 2.1. \21l \2‘2!j fweak divergence) The weak divergence of v G V{K), 
denoted by • v, is a hounded linear functional in the Sobolev space so that 

its action on any (f G H^{K) is given by 

(Vu, • V, (j))K ■■= -(vo, V<j))K + (Vb • n, (l))dK, 

where n is the unit outward normal direction on dK. 

For any non-negative integer r, denote by Pr{K) the set of polynomials with 
degree r or less on K. 

Definition 2.2. f21[ \22l/ f discrete weak divergence) The discrete weak divergence 
of w G V{K), denoted by '^w,r,K ■ v, is the unique polynomial in Pr{K), satisfying 

(2.1) • V, (/))k =-(vo, -b (vb • n, \/(j)GPr{K), 

where n is the unit outward normal direction on dK. 

Definition 2.3. WR\^ fweak gradient) The weak gradient ofv G V{K), denoted 
by VujV, is a matrix-valued bounded linear functional in the Sobolev space [H^(K)]‘^^‘^, 
so that its action on any (p G [H^{K)]‘^^‘^ is given by 

(Vu,v, p)k ■■= -(vo, V • (p)k + (vb, pn)gK, 

where n is the unit outward normal direction on dK. Here the divergence SJ ■ p is 
applied to each row of p, and pn. is the usual matrix-vector multiplication. 

Definition 2.4. \2^ (discrete weak gradient) The discrete weak gradient of 

V G V{K), denoted by Vw,r,K'v, is the unique matrix-valued polynomial in [Pr{K)]‘^^‘^, 
satisfying 

(2.2) = -{^0,^ ■ t)k + {Vb,Tri)9K, \/(p G [PriK)f^‘^, 
where n is the unit outward normal direction on dK. 
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Using the weak gradient, we may define the weak strain tensor as follows 

(2.3) eu,(u) = i(Vu,u + Vu,u^). 

Analogously, the weak stress tensor can be defined by 

(2.4) cru.(u) = 2/x£u,(u) + A(Vu, • u)I. 

3. Numerical Algorithms. Let 7)i be a finite element partition of the domain 
U C consisting of polygons in 2D or polyhedra in 3D which are shape regular 
interpreted as in [21]. For each T G Th, denote by hx the diameter of T. The mesh 
size of 7h is defined as h = ma^TeTh Let T G 7/t be an element with e as an edge 
in 2D or a face in 3D. Denote by £h the set of all edges or faces in Th and = £h\dfl 
the set of all interior edges or faces in Th- 

On each element T gTh, denote by RM(T) the space of rigid motions on T given 
by 


RM(T) = {a + ?7x : a G r] G so{d)}, 

where x is the position vector on T and so{d) is the space of skew-symmetric d x d 
matrices. The trace of the rigid motion on each edge e C T forms a finite dimensional 
space denoted by 

Prm(.g) = {v G [L^{e)Y : v = v|e for some v G RM(T), e C dT}. 

For any positive integer fc > 1 and element T £ Th, we introduce a local weak 
finite element space as follows 

(3.1) U(fc,T) = {v = {vo,v4 : Vq G [Pfe(T)]^ v, G Ufe_i(e)} , 

where Vfc_i(e) = [Pfe_i(e)]‘^ + PRM(e). Since PRuie) C Pi{e), then the boundary 
component 14_i(e) is given by [Pk-i{e)Y for fc > 1 and PRM{e) for k = 1. 

The global weak finite element space 14, is given by patching the local spaces 
V{k,T) with a single-valued component v;, on each interior element interface. All the 
weak finite element functions v G 14 with vanishing boundary value v;, = 0 on F form 
a subspace of 14, which is denoted as 

(3.2) 14° = {v = {vojVb} G 14 : Vf, = 0 on F}. 

For each v G 14, the discrete weak divergence Vw^k-i ■ v and the discrete weak 
gradient Vu;,fc-iv are computed by using (EH) and (12.21) on each element T G 74; i.e., 

(Vu,,fe-1 • v)|7’ =V.u],k-l,T • (v|t), V G 14, 

(Vu,,fe-iv)|7’ =Vu,^fe_i^T(v|T), V G 14. 

For simplicity of notation, we shall drop the subscript k—1 from the notation 
and in the rest of the paper. For each edge or face e G £h, denote by Qh the 

projection operator onto the space 14_i(e); i.e., the polynomial space [Pfe_i(e)]^ 
for fc > 1 or the rigid motion space PrmY) for fc = 1. 
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Next, we introduce two bilinear forms 
(3.3) s(w,v)= ^ h:^^{QbWo-MVb,QbVo-Vb)dT 

T€Th 


(3.4) as(w,v)= ^ 2(/x£i„(w),eii,(v))T + ^ (AV„ ■ w, • v)t + s(w, v), 


TGTfc 


TdTh 


where 6^ (w) and V^, • w are computed by using the discrete weak gradient and weak 
divergence operators. 

Weak Galerkin Algorithm 1. For a numerical solution of the elasticity 
problem U.3\) . find u/j = {ug, Uf,} £ Vh with u;, = on T such that 


as(u/,,v) = (f,Vo), Vv = {vo,v4 € V^°. 


(3.5) 


The rest of this section is devoted to a study of (13.51) on the solution existence 
and uniqueness. 

Theorem 3.1. There exists one and only one solution to the weak Galerkin finite 
element scheme (TOP . 

Proof. Since the number of equations equals the number of unknowns in (13.51) , it 
suffices to prove the solution uniqueness. To this end, let G 14, j = 

1, 2, be two solutions of (13.51) . It follows that = (5f,u on T and 


= (fWo), Vv = {vo,Vf,} G V)°, 3 = 1,2. 

The difference of the two solutions, 'w — satisfies w and 


(3.6) 


as(w,v)=0, Vv = {vo,v4 G 14°. 


By letting v = -w in (13.61) we obtain 


as(w, w) = 0. 


Using the definition of as(-, •) we arrive at 



Ten 


Ten 



Ten 


which implies that 


(3.7) 

(3.8) 


w(w)=0, in T, 

(5t,wo - Wf, = 0, on dT. 


From the definition of the weak gradient, we have 

(V„w, t)t = (Vv^fQ, t)t - (wq - Wb, Tn)OT 

= (V-wo, t)t - {Qb'ffo - 'Wb, rn)aT 
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for all T S [Pk-i{T)Y^'^,k > 1. It follows from (13.81) that Vwq = Vu,w on each 
element T. Thus, with the help of (13.7L 

e(wo) = eu,(w) = 0, 

which leads to wq G RM(T) C [Pi{T)]^. It follows that Wo|e = Q&wq = Wf,, and 
hence wq is a continuous function in with vanishing boundary value on F. From 
the second Korn’s inequality (19.1311 . we obtain wq = 0 in II, and hence Wb = 0 from 
(13.81) . This shows that and hence the solution uniqueness and existence. 

We remark that the result holds true for any A > 0. □ 

4. An Equivalent Mixed Formulation. A strong form of the mixed formula¬ 
tion (ll.4ll - (ll.5ll reads as follows: Find u and p satisfying u = u on F, the compatibility 
condition X~^pdx = /p u • nds, and the following generalized Stokes equations 

^ -V • (2^£(u)) -P Vp = f, in II, 

V • u = X~^p, in H. 

We assume that the generalized Stokes problem (I4.1|l has the i7^+®(II) x i7®(II)- 
regularity for some s G (^, 1] in the sense that, for smooth data f and u, the solution 

u and p of (El) satisfy u G [77^+'’(II)]^^, p G i7®(II), and the following a priori estimate 

(4-2) l|u||i+s -f IIpIU < C'(l|f|U-i + l|uL+i,r) 

for some constant C independent of the parameter A. The regularity estimate (14.21) 
can be found in [IQIIT] for convex polygonal domain with s = 1 when A = oo. For 
large values of A, one may heuristically apply the result for the Stokes (by viewing 
A“^p as a given function) to obtain 

\Mi+s + IblU < C'dIflU-i + l|ub+i,r + A”dblU), 

which implies (14.2p when A is sufficiently large. 

The idea of weak Galerkin can be applied to the mixed formulation ()1.4l) - (ll.5p 
for the linear elasticity problem (ICT-lfO]). This application requires an additional 
finite element space that approximates the auxiliary variable p. More precisely, we 
introduce 


Wh = {q: q\T G Pfc-i(r), T G %}, = Wh H PHn) 

and the following bilinear forms 

a(’w, v) = 2(p£u;(w), ew(y))h + s(w, v), 

6(v,g) = (Vu, ■v,q)h, 
d{p,q) = X~^{p,q), 

where w, v G Vh, p,q & Wb, and 

(p£u,(w),eu,(v))b = ^ (Ai£»(w),£u,(v))T, 

TdTh 

(Vu, ■ V, q)h = ^ (V„ • V, q)T. 

TdTh 
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Weak Galerkin Algorithm 2. For a numerical solution of the linear elasticity 
problem ^1.4\ )- U.5\) . find U/j = {uq, U;,} G I 4 and p^ G Wh satisfying Uf, = on F, 
the compatibility condition {X~^pfi, 1 ) = /p u • nds, and the following equations 

(4.3) a(u^,v) + 6(v,p?i) = (f,Vo), Vv G 

(4.4) b{\ih, q) - d{ph, q) = 0, Vg G 


Lemma 4.1. The weak Galerkin Algorithms'^ and\^ are equivalent in the sense 
that the solution U/j from \3.5\) and i4-3p - (T^ are identical to each other. 

Proof Assume that u/i and pu solves (I4.3D - (I4.4I) . Note that the equation (14.41) 
can be rewritten as 

(V^ ■ Uh,q)T - X~^{ph,q)T =0, Vg G Pk-i{T). 

Since • u/j G Pk-i{T), then pu can be solved from the above equation as 
(4.5) ph = X\7.u] ■ Uh- 

Substituting (14.51) into (14.31) yields 

a(u/,,v) + A6(v, V„ • Uh) = (f,vo), Vv G V°. 


Note that 


a(w,v) + A 6 (v, •-w) = as(w,v), Vw, v G 14 . 

Thus, we have 

as(u^,v) = (f,vo), Vv G 14°, 

which is the same as (13.51) . It then follows from the solution uniqueness that Uh is 
identical with the numerical solution arising from the weak Galerkin Algorithm [ 1 ] 

A similar argument can be applied to show that if solves (13.5L then the pair 
(u/i; XVw • u?i) is a solution of (I4.3I) - (I4.4I) . Details are left to interested readers as an 
exercise. 0 

The reformulation (fT:i)-(IT31) of the elasticity problem as a generalized Stokes 
system with nonzero divergence constraint often leads to numerical approximations 
which are locking-free as A —> 00 . Due to the equivalence between the primal formula¬ 
tion (13.51) and the mixed formulation (14.31) - (14.41) in the WG finite element method, the 
WG scheme in the primal formulation (13.51) is locking-free if the corresponding WG h- 
nite element method for dEl-dlSl) can be proved to be stable and accurate in terms of 
the parameter A. The rest of the paper is devoted to a stability analysis for the mixed 
weak Galerkin finite element scheme ()4.3D - (I4.4I) . which in turn implies a locking-free 
convergence for the linear elasticity problem in the displacement formulation (13.51) . 

5. Stability Conditions. In the weak finite element space 14, we introduce the 
following semi-norm 

v||| = ( l|£(w)|lT + /i?iQf,Vo-Vb|||p.y, 

TdTh 


(5.1) 


V G 14 . 
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Lemma 5.1. The semi-norm ||| • |||, as defined by is truly a norm in the 

linear space V^. 

Proof. We shall only verify the positivity property for ||| • |||. To this end, assume 
that III V III = 0 for some v € Vf^. It follows that e(vo) = 0 on each element T G Th and 
Qb'Vo = Vb on dT. Thus, Vq £ RM(T) satisfies Vo|e = Qb(vo|e) = v^, which implies 
the continuity of Vq in the whole domain 12. The boundary condition Vf, = 0 on T 
implies that Vq = 0 on T. From the second Korn’s inequality (I9.13L we obtain vq = 0 
in 12. The fact that Qb'Vo = vt on dT gives vt = 0 on dT. Thus, v = 0 in 12, which 
completes the proof of the lemma. □ 

Lemma 5.2. There exist positive constants ai and 02 such that 
(5.2) ai|||v|f < a(v,v) < a 2 |||v|f, V v G I4°- 


Proof. From (12.21) and the integration by parts, we have 

(e„(v),(/?)T =(Vvo,i(v3 + V5^))T - (vo -Vb,h(f-\-ip'^)n)9T 
( 5 . 3 ) 2 ^ 2 

= (e^(vo), V?)t - {Qb^o - Vb, 2 + V5^)n)aT 

for any ip G [Pk-iiT)]‘^^‘^. Thus, 

(5.4) |(e„(v),</j)T| < ||e(vo)||T||i^||T + HQbVo - Vb||aT||^(‘/? + <p'^)n||aT- 
From the trace inequality (19.21) and the usual inverse inequality we have 

\\^{ip + ip'^)n\\9T < C'(/i;^^||(p||^ + /it||V(^||^)" 

<chpMT. 

Substituting the above into (15.41) yields 

|(e.^(v),(/?)T| < (||e(vo)||r + C'h^"||QbVo - vbllar) \W\\t: 
which leads to 

(5.5) lkj«(v)||T < 2 ||£(vo)||| + ChffWQb^ro - Vb||ar. 

By representing (e(vo),i^)T in terms of other two terms in (15.3|) . we can derive 
the following analogy of (15.5() 

(5.6) lk(vo)||T < 2||e„(v)||y + Ch'f.^WQb'VQ - ^rb\\^T■ 

The left inequality in (j5.2() is a result of (15.6p by summing over all the element T gTh, 
and the right inequality can be obtained by summing (15.5p over T G T/t- D 

In the finite element space we introduce the following norm 

Ml = \<l\lb + h^yq\\lb, 


qeWh, 
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where 


\<l\lH=hY.\Me\\l \\yq\\lH= E ll^'^ll 

eGS2 TGTh 


2 

T- 


(5.7) 


Lemma 5.3. There exists a constant /3 > 0 such that 

b{yr,q) 


sup 


> /?|||g|||o, V q e 


Proof. From the definition of the discrete weak divergence (EH), we have 

(5.8) 6(v, g) = E = E {“(’^ 0 , Vg)T + (v^ • n, q)oT} ■ 

Ten Ten 

By setting v = v, := {—/i|g]en} in (15.81) . where e G f° and n is the unit 
outward normal direction on e, we arrive at 

(5.9) 6(v„9) = h^Y. I|V9|It + /^ E IIMelle = llklllo- 

Ten eeeo 

Furthermore, it is not hard to see that there exists a constant Cq such that 


(5.10) 

Thus, we have 


sup 

vGVft°.V5^0 


ll|v,||| < Collklllo. 
H^,q) ^ bf^q^q) _ Ilk 


> 


>t^o-'ll 


IIO’ 


which proves the inf-sup condition (I5J1) . □ 


6. Preparation for Error Estimates. For each element T G 7)i, let Qo be the 
projection onto [Pk{T)]'^. For each edge/face e C dT, recall that Qt is the 
projection onto the finite element space on e. Denote by Qh^i the projection onto 
the weak finite element space 14 such that on each element T € Th, 

QhU := {Qou, Qfeu}. 

Furthermore, let Qh and be the projection onto Pk-i(T) and [Pfc_i(T)]'^^'^, 
respectively. 

Lemma 6.1. JMEI For any v G [77^(D)]‘^, the following identities hold true for the 
projeetion operators Qh, Qh, o,T,d Qh 


(6.1) V^-(Q,,v)=Q^(V-v), 

(6.2) V^(Q„v) =Q;,(Vv). 


Consequently, one has 
(6.3) 


SwiQh'^) = Qh^i^)- 






11 


Proof. We outline a proof for (16.11) only; a similar approach can be adopted to 
prove (16.21) . From (12.11) and the integration by parts we have 


(V^ • (Q/tv), (^)t = -(Qov, V(^)t + (QbV • n, if)oT 
= -(v, V(^)t + (v • n, if)oT 
= (V • 

= (Q?i(V ■ 

for all ip e Pk-i(T). Since by construction • ((5?iV) e Pk-i(T), then (16.11) follows. 

□ 


Lemma 6.2. Assume that (w;p) € x 7 > satisfies the 

following equation 


(6.4) 2V • (/ie(w)) + Vp = — 77 , in fl. 

Let [Qh'W] Qhp) be the projection of (w; p) in the finite element space 14 x 114 • 
Then, we have 

(6.5) 2{peni{Qh'v^), £w{'v))h + (Vu, • v, Qhp)h = {v, w) + ^w(v) + 6»p(v), 

for all v G , where £w and Op are two functionals in the linear space given by 


(6.6) £w(v) = 2 ^ (vo - Vb,p(e(vif) - Q^e(w))n)aT, 

TdTh 

(6.7) 6»p(v) = X! “ Q?iP)n)aT- 

TdTh 


Proof From (16.31) . (12.21) and the integration by parts, we get 


2p{e^,{Qhy^),£wM)T 

=2p(Q?,£(w),e^(v))T 

= - 2p(vo, V • (Q^e:(w)))T + 2p(vb, Clf^e{w)n)oT 
=2p(Vvo, Qb£(w))T - 2p(vo - Vb, Qb£(w)n)aT 
=2p(Vvo, £(w))t - 2p(vo - Vb, Qb£(w)n)aT- 


( 6 . 8 ) 
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Next, we have from (EH), the integration by parts, and P^)dT = 0 that 

(V^ • V, Qhp)h = (V„ • V, Qhp)T 

TGTh 

= {-i^o,^iQhp))T + {vbAQhp)n)aT} 
tgU 

= ^ {(V • Vo, Qhp)T - (vo - Vb, {Qhp)n.)oT} 

T&Th 

= X! {(^ • '"0, p)t - (vo - Vb, (Qbp)n)aT} 

TgTh 

= X! {“('"0, Vp)t + (vo, pn)oT - (vo - Vb, {Qhp)n)oT} 

T&n 

= ^ {-(vo, Vp)t + (vo - Vb, pn)oT - (vo - Vb, iQhp)n)dT} 

TGTh 

= - (vo, Vp) + X! “ Qhp)n)oT, 

Ten 

which leads to 

(6.9) (vo, Vp) = -(Vu, • V, Qhp)h + ^ {'^0- Vb, (p - Qhp)n)oT. 

TdTh 

Now testing (16.41) by using the component vq of v = {vq, Vb} € V))* yields 

(6.10) - 2(V ■ (p£(vif)), Vo) - (Vp, Vo) = (rj, vq). 

From the integration by parts, we can rewrite (I6.10p as 

(6.11) 2 '^vo)t - 2 X! (/^eWn, vo)aT - (Vp, vo) = (r?, vq). 

TdTh TGTh 

Substituting (16.81) and (16.91) into (16.111) yields 

2p X! + (vo - Vb,Qbe(w)n)aT - (vo,e(w)n)aT} 

T&Th 

+ (Vu, • V, Qbp)/i - X! ^'"0 “ ~ Qhp)ri)dT = {rj, vq), 

TeTfe 

which implies that 

2p X] £w{-v))t - (vo - Vb, (£(w) - Qb£(w))n)aT - (vb, e(w)n)aT} 

TgTh 

+ (Vu, • V, Qhp)h - {^0- Vb, (p - Qhp)n)oT = {r], vq). 

TGTh 

Using the boundary condition Vb = 0 we obtain 

2(pe„(Qbw),e.u,(v))b + (V^ • v, Qhp)h 

(6.12) = M(e(w) - Qbe(w))n)aT + (vo - Vb, (p - Qhp)n)oT} 

TGTh 

+ (^7,vo), 

which is precisely the equation (16.5F This completes the proof. □ 
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7. Error Estimate in a Discrete i7^-Norm. For the weak Galerkin finite 
element solution {\ih]Ph) = ({uq,G 14 x Wh arising from (14.31) - (14.41) for the 
linear elasticity problem (fTH-ifr^. we define its error functions Gh and 4 by 

(7.1) Gh = {eo, et,} = {Qou - uq, - Ub}, 

(7.2) Ch = QhP-Ph, 

where (u;p) is the exact solution of the variational problem (ll.4l) - (ll.5l) . It is clear 
that Gh € V° and C/t G W°. 

Lemma 7.1. The error functions Gh and fh defined in satisfy the 

following error equations 

(7.3) a(e;,,v) -h &(v,4) = V?u,p(v), Vv G 14°, 

(7.4) b{Gh,q) - d{Ch,q) =0, Vg G 

where 


PM,piy) = ■^u(v) -h 6(p(v) -h s(QbU, v). 


Proof. Observe that the exact solution (u;p) satisfies the equation (16.41) with 
rj = i. Thus, from Lemma [6. 2 1 we have 

2{p£nj{QhU),£ni{v))h + (Vu, • V, Qhp)h = (f, Vq) -|- ^u(v) -|- 6»p(v), 

which leads to 

(7.5) a{QhU, v) 6 (v, Qhp) = (f, vq) ^u(v) -h 6»p(v) s{QhU, v). 

Subtracting (1431) from dig gives the equation (1731) . 

To derive (17.4L using (16.11) we have for any q G 114 

g (Vu, • iQhu),q) - X~^iQhP,q) = (Qhiy ■ u),g) - X~\QhP.q) 

= (y ■u,q)- X~^{p,q) = 0, 

where we have used dESl) in the second line. The difference of (17.61) and (14.41) yields 

(1731) . □ 

We are now in a position to derive an error estimate for the weak Galerkin finite 
element approximation {uh]Ph)- 

Theorem 7.2. Let the solution of o-ra be sufficiently smooth such that 
(u;p) G [77^+^ (O)]"^ X 77^(0) for some k > 1. For the weak Galerkin finite element 
solution {uh;ph) G 14 x 114 arising from we have 

(7.7) |||Q/«u-u/,||| -1 A“5||g^p-p^|| + \\\Qhp-ph\\\o < Cli''(||u||fc+i ||p||fc), 

where C is a generic constant independent o/(u;p). Consequently, the following error 
estimate holds true 

u - Uh\\\ + x~i Up -Ph\\ + lb - Phio < Ch'^{\\u\\k+i + IblU). 


(7.8) 
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Proof. By choosing v = e/j in (17.31) and g = in (17.41) we have 
n(e/i, e/i) A ||C/i|| “ • 

By applying Lemma |9.11l to the term (/?u,p(e?i) we arrive at 
a(e;„e^) + A-i||a|P <C'/r'=(||u|| 

fe+i + lblU)ll|e.|||. 

Next, using Lemma 15.21 and the above estimate we obtain 

(7.9) a,|||e,f + A-i||af < Ch>^{\\nh+i + |blU)ll|e/.|||. 

To derive an error estimate for the “pressure” variable p in a A-independent norm, 
we use the inf-sup condition (15.71) to obtain 

(7.10) /llliaillo < sup 

vev“,v5^o llrll 

From dlSl) 


&(v, Ch) = -a{eh, v) + Pu,p(v). 

Thus, it follows from Lemma [Ql the error estimate (EH), and Lemma [9.111 that 

|&(v,C?.)| < a2|||e,,||| |||v||| + |pu.p(v)| 

< Ch'=(||u||fc+i + IblU) |||v|||. 

Substituting the above estimate into (I7.10|) yields 

(7.11) iiiaiiio<c'^Miuiu+i + iHU). 

Combining (17.91) with (I7.11|) gives rise to the error estimate (17.71) . Finally, ()7.8p stems 
from the usual triangle inequality, the estimate (TTtI) and the error estimate for 
projections. □ 

8. Error Estimate in As usual, we use the duality argument to derive an 
error estimate for the weak Galerkin finite element method. To this end, consider 
the problem of seeking $ € [iL^(fl)]'^ and ^ G ^^(n) satisfying 


(8.1) 

-V.(2/r£($)+e/)=eo, 

in fl. 

(8.2) 

V.$-A-^^ = 0, 

in n. 

(8.3) 

$ = 0, 

on F. 


Assume that the dual problem (I8.1D - (I8.3I) has the x i7®(r2)-regularity with 

A < s < 1 in the sense that the solution (^;'C) £ [iL^"'"®(S4)]^ x iL®(S4) and satishes 
the following a priori estimate: 

(8.4) ||$|U+i + ||^|U<C||eo||. 

Theorem 8.1. Assume that the solutions of (EIF are suffieiently smooth such 
that (u;p) G x for some integer k > 1. Let {uh;ph) G 14 x Wh 
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be the corresponding weak Galerkin finite element solution arising from !i4.3\ )- i4.4\ )- 
Then, under the regularity assumption iS.4\l , there exists a constant C, such that 

(8.5) IIQou-uoll < C'/i''+®(||u||fe+i + ||p||fc). 

Moreover, it follows from the triangle inequality and the error estimate for the 
projection that 

(8.6) ||u - uoll < Ch’^+‘’ {\\u\\k+i + IblU). 


Proof. Note that the solution (^;^) of (I8.1|) - (I8.3I) satisfies (16. 4|) with rj = eg. 
Thus, using Lemma [121 namely the identity (16.51) . we have 

(8.7) 2{iJ,e.u]{Qh^),£w{'v))h + (Vu, • V, Qhf)h = (eg, vq) +^$(v) + %(v) 
for all V e 

By choosing v = e^, in (18.71) we obtain 

( 8 . 8 ) ||eg|p = a{Qh^, e^) + b{eh, Qhf) - 
where 

Using dUD we have 

(8.9) &(e/^,Q;^0-'^(2/^e,C/^)=0. 

From (Hr]) and (|H21) 

b(Q/,^,a)= ^(V„-Q,,$,a)T 

TGTh 

= ^(Q/.(V-$),a)T 

Ten 

= d{QU,a). 

Thus, by substituting the above into (18.91) 

b{eh,Qhf) = b{Qh^,Ch). 

Combining the last equation with (18.81) we obtain 

I|egf = a{Qh^,eh) + b{Qh^,Ch) - 
which, together with the error equation (17.3L leads to 

(8.10) ||eo|p = (fu,p{Qh^) - 

The rest of the proof shall deal with the two terms on the right-hand side of 
(I8.10|) . The second term (pu^q{eh) can be estimated by using Lemma [9.111 the error 
estimate (1721) . and the regularity assumption (18.41) as follows 

l¥>$,j(e^)| ^Ch^dl^lU+i + lieiU) IKIII 

<C'h'=+^(||u|U+i + |b|U) ||eo||. 


( 8 . 11 ) 






16 


The estimate for the first term <fu,piQh^) is ^ bit complicated, for which we detail 
as follows. 

(i) For the term s((5/iU, Qt^) in y^u.piQh^), from the Cauchy-Schwarz inequality, 
(19.11) and (19.31) . we have 

|s(Q;iU, =1 ^ {QbiQou) - QbU,Qb{QQ^) - Qb^)dT 

TdTh 

TeTh TGTfe 

<C/r'=+^||u|U+i ms+i- 


( 8 . 12 ) 


< 




(ii) For the term £u{Qh^), we have from the orthogonality of Qb and the boundary 
condition (18.3|) that 

(e(u) - Q^^(e(u)))n)^T 

Ten 

=2fi ^ {Qb^ - $,e(u)n)aT = 0. 

Ten 

Thus, it follows from the Cauchy-Schwarz inequality, dUD, (USD and (19.4|1 that 


(8.13) 


\£^{Qh^)\ = 




dT 


Ten 


(8.14) 


2m XI ^*30^ - (e(u) - Q^e(u))n) 


dT 


Ten 




dT I 


<2m( X ^T||e(u) - Q;je(u)|||j,) ' ( X ^t^IIQc 
Ten Ten 

<Ch^+^\\u\\k+i ||$|U+1. 


(iii) As to the term 9p{Qh^), we again use the orthogonality of Qb and the 
boundary condition (18.31) combined with the Cauchy-Schwarz inequality, (19.11) . (19.31) 
and (ITO to obtain 


\0p{Qh^)\ = 


X {Qo^-Qb^,{p-QhP)r^) 


dT 


Ten 


(8.15) 


X (Qo^ - (P - QhP)n) 


dT 


Ten 


< 


0^- ^lloT 


( X hrWp-QhPWlrW X ^t^WQ' 

Ten Ten 

<ch^+^\\p\\k ||$iu+i. 

By combining the three estimates (18.121) . (I8.14p . and (I8.15P we arrive at 

I<Mu.p(Q/.$)| < Ch’^+m\uU+i + lblU)||^IU+i 

<C/r'=+^(||u|U+i + |b|U)||eo||, 


(8.16) 
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where we have used the regularity assumption (18.41) in the second line. Finally, by 
substituting (18.111) and (18.161) into (18.101) we obtain the desired error estimate of (18.51) . 
This completes the derivation of the error estimate. □ 

9. Supporting Tools and Inequalities. In this section, we present some tech¬ 
nical inequalities that support the error analysis established in previous sections. 

Recall that 7^ is a shape-regular finite element partition of fl. There exists a 
constant C > 0 such that, for any side e C T with T € Th, the following trace 
inequality holds true m 

(9.1) Ml<c{h^^ufT + hTmfT), yct>eH\T), 

where /it is the size of T. Furthermore, in the polynomial space Pj{T), j > 0, we 
have from the inverse inequality that 

(9.2) ||<^||2 <c/i^i||0ll|, v0GP,(r), Ter^. 

Lemma 9.1. Assume that 7h is a finite element partition ofD, satisfying the 
shape regularity assumption as defined in \2Vf . Let k > 1 be the order of the finite 
element method and 1 < r < k. Let w G [77’'“*'^(0)]'^, p G H^{n) and 0 < m < 1. 
There holds 

(9.3) ^ /.?P||w-gow||^,„<C'/i^(’-+i)||w||2+i, 

TeTh 

(9.4) ^T"lk(w) - Qh£i^)\\T,m < C'/l^1|w||2_^i, 

Ten 

(9.5) ^ h^T^\\p-QHprT^^<Ch^^\\p\\l. 

TeTh 

9.1. Korn’s inequality. Korn’s inequality is a fundamental tool in the study 
of elasticity equations. The inequality can be found in many existing literature, see 
[mill [mil] for example. For convenience, we provide a summary here for this useful 
inequality. 

Theorem 9.2. (Korn’s Inequality) Assume that the domain LL is open bounded 
with Lipschitz continuous boundary. Then, there exists a constant C such that 

(9.6) ||Vv||o<C(||e(v)||o + ||v||o), V v G [^'(fl)^ 


Proof. A proof can be given by using the following inequality m 

(9.7) ||<^||o<C'(||V<^||_i + ||<^||_i), V</.gL2(L!). 

To this end, for any v G [77^(11)]'^, it is not hard to check that 

(9.8) djdkVi = dj£ik{-v) + dk£ij{v) - diEjki^). 

It follows that Vui G if e(v) G Moreover, from (19.71) and (19.8|) 

one has 


||Vv||o<C(||V£(v)||_i + ||Vv||_i) 
< C'(lk(v)llo + ||v||o), 








18 


which is the Korn’s inequality ( 1 ^ . □ 

The following is a characterization of the space of rigid motions as the kernel of 
the strain tensor operator. 

Lemma 9.3. Let Lt be an open bounded and eonnected domain in M'’*. The space 
of rigid motions RM(f2) is identical to the kernel of the strain tensor operator; i.e., 
for V £ [H^(Ll)]'^, there holds e(v) = 0 if and only ifv€ RM(n). 

Proof. For any v £ RM(n), there exist a £ and a skew-symmetric dxd matrix 
77 such that u = a -I- px. It is easy to check that e(v) = 0. 

For any v satisfying £(v) = 0, we have from (19.811 that djdkVi = 0, and hence 
V £ [Pi(n)]'^. Thus, there exist a £ and 77 £ such that v = a-|- 77X. Since 

e(v) = 0 , we get 

0 = £(v) = £(a -I- rjx) = ■i(V(a -f rjx) + V(a -f 77 x)’^) = ^(77 -|- 77 '^). 

It follows that 77 ^ = — 77 , which means that 77 is skew-symmetric, and hence v £ RM(0). 

□ 

Lemma 9.4. Let x^, i = 0, ...,d — 1, be d-points in LI which form a {d — 
1) —dimensional hyperplane. If 'v G RM(n) satisfies v(xi) = 0, i = 0,..., d — I, then 
we must have v = 0. 

Proof. Let v = a -f 77X be the rigid motion with v(xi) = 0. Thus, 

(9.9) 77(xi-xo) = 0, i = l,...,d-l. 

Note that the set of vectors {x^ — XQ}f~l is linearly independent since they form a 
basis of a (d — I)—dimensional subspace of R'^. 

For d = 2, the skew-symmetric matrix 77 is either zero or invertible. From the 
equation (19.911 . we see that 77 can be nothing except 77 = 0. For d = 3, the matrix 77 
has eigenvalue Aq = 0. If 77 7 ^ 0, then the eigen-space corresponding to the eigenvalue 
Aq = 0 must have dimension 1. But the equation (19.91) indicates that the dimension 
for this eigen-space is no less than 2. Consequently, we must have 77 = 0. Finally, 
from v(xo) = 0 we have a = 0. This shows that v = 0. □ 

Theorem 9.5. (Second Korn’s Inequality) Assume that the domain LI is con¬ 
nected, open bounded with Lipschitz continuous boundary. Let 4> : H^{L1) —> R+ be a 
semi-norm on H^{L1) satisfying 

<I>(v) = 0 and v £ RM(n) => v = 0. 

Then, there exists a constant C such that 

(9.10) ||v||i < C(||£(v)||o+ $(v)), 
for alive [-ff^(f^)]'^- 

Proof. We verify the inequality (19.101) by a contradiction argument. To this 
end, assume that (19.101) does not hold true. Then, for each integer n, there exists 
v„ £ [H^{Ll)]‘^ such that 


( 9 . 11 ) 


V„||l > 7r(||£(v„)||o -f ^(Vn)). 
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We may assume ||v„||i = 1, and hence, there exists a subsequence {v„j,} which is 
weakly convergent in and strongly convergent in From (19.111) . we have 

||e(vn)||o + ^(v„) < n~^. 


Thus, 


lk(Vn)||o + $(V„) 0, 

which, together with Korn’s inequality (19.6L implies that {v^^,} is a Cauchy sequence 
in Hence, there exists a function v S such that 

Vrifc —>■ V strongly in 


Moreover, we have 


||e(v)||o + $(v) = lim (||e(v„J||o + $(v„J) = 0. 

K—^OO 


Thus, 


V G RM(H) and ^(v) = 0, 

which leads to v = 0. This is a contradiction to the assumption that 

||v||i = lim ||v„J|i = 1. 


This completes the proof. □ 

The following are two particular cases of the seminorm $(•) that satisfy the con¬ 
ditions of Theorem 19.51 

Corollary 9.6. Assume that the domain H is connected, open bounded with 
Lipschitz continuous boundary. Let D C LI be a subdomain of LI, and Qd.rm be 
the projection from [L^{L1)]‘^ onto the space of rigid motions RM(Zl). Then, there 
exists a constant C such that 

(9.12) ||v||i < C'(||e:(v)||o + \\Qd,rm'v\\o,d), 

for all V G 

Proof. Define 5*(v) = ||(5 _d,rmv||o,d which is clearly a semi-norm in 
The estimate (|9.12l) stems from Theorem 19.51 if the conditions of the theorem are 
verified for this semi-norm. To this end, for any v satisfying £(v) = 0, we have from 
Lemma [9.31 that v G RM(D), and hence Qd.rm'v = v a.s D C LI. If, in addition, 
^(v) = 0, then ||v||o,d = \\Qd,rm^\\o,d = 0. Hence, v = 0 in H. □ 

Corollary 9.7. Assume that the domain LI is connected, open bounded with 
Lipschitz continuous boundary. Let Fi C dLl be a nontrivial portion of the boundary 
dLl with dimension d—\. For any fixed real number 1 < p < oo, there exists a constant 
C such that 


(9.13) 

for all V G [i7^(D)]‘^. 


||v||i < C(||e(v)||o -f ||v||ip(ri)). 
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Proof. In we define a semi-norm by 

‘^(v) := (^ |vps) ", V e 

We claim that this semi-norm satisfies the conditions of Theorem In fact, if V 
satisfies 5*(v) = 0, then we have v = 0 on Ti. Hence v = 0 on d-points of a hyperplane 
of dimension d — 1. If, in addition, v £ RM(H), then from Lemma 19.41 we must have 
V = 0 on H. □ 

9.2. Some technical inequalities. For any xq G denote by i?(xo,r) the 
d-ball centered at xq with radius r. The unit d-ball centered at the origin is denoted 
by B, called the reference d-ball. The reference d-ball can be identified with the d-ball 
i3(xo,r) using the following affine map 

X = F(x) := Xq -I- rx. 

The inverse of the affine map F is given by 

X = F"\x) := (x - Xo)/-r. 

Any function (p = (pfx.) on B{xQ,r) defines a function on the reference d-ball as follows 

p{x) = (/)(F(x)), 

which shall be denoted as 

<^(x) = </)(x). 

It is clear that </> £ RM(i?(xo, r)) if and only if ^ £ RM(i3). 

Lemma 9.8. Let QB{xo,r) Qg projections onto the space of rigid 

motions RM(i3(xo,r)) and RM(i3), respectively. Then, for any v £ [L^(i?(xo,r))]'^ 
we have 

(9.14) (3B(xo.r)V = Qgv- 


Proof. From the definition of the projection, we have for any (p £ RM(i3(xo, r)) 


/B(xo,r-) 


QB(xo,r)^ P^dx 


IB(xo,r) 


V pdx = / V pdx = r‘^ Qb^ 

Jb Jb 


pdx. 


Also by changing the domain from B{xQ,r) to B, 

/ QB(xo,r)V (pdx = r‘^ I QB(xo,r)V pdx. 

J B{fx.Q,r) JB 


fB 


Q B{x.Q,r)^ / Q(f)dx^ 


It follows that 
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which leads to (|9.14l) . □ 

Lemma 9.9. There exists a constant C such that, for any w S [H^{B{x.Q,r))]‘^ , 
we have 

||w Qi?(xo,r) w||o,B(xo,r) ^ || 0,i?(xo ,r) • 


Proof. Let = w — QB(xo.r)W. It follows from Lemma [9.81 that 

= w — QjjW. 

It is also easy to see the following identities 

(9.15) = 0, e(w^) = £r(w). 

By mapping to the reference d-ball, we have 

llw^llo.B(xo.r) 

Using the second Korn’s inequality (19.121) with ft = B and D = B, we obtain 


r-L||2 




< C\\e{w)\\l g 

where we have used (I9.15|) in the second line. Thus, one has 

llo.B(xo.r-) < lk(w) ||o,B(xo.r-) • 

This completes the proof of the lemma. □ 

Lemma 9.10. Let It be a shape regular finite element partition ofQ. There exists 
a constant C independent of T G Th such that 

(9.16) ||vo - QbMlr < ChT\\e{yo)\\l, V vq G [PkiT)f. 


Proof. From the shape regularity assumption, there exists a d-ball i?(xo,r) C T 
for which the radius r is proportional to Lt] i.e., r = Aq/it with a constant Aq 
bounded away from 0. For any Vq G [PkiT)Y, consider the projection QB(xo.r)Vo 
which is naturally extended to T. Since {QB{xo,r)'Vo)\dT belongs to the finite element 
space on dT and Qb is the projection onto this finite element space, then 

(9.17) ||vo - Qb^oWar < ||vo - (3B(xo.r)Vo||aT- 

Using the trace and inverse inequality ()9.1I1 - (I9.2|) we have 

||vo - Qb{ Xq, r)Vo||aT < C{hrp^\\vo — QB(xo,r)Vo||T + /it||V(vo — Q_B(xo.r-)Vo) ||t) 

< Ch-^\\vo - (3B(xo,r)Vo||T 

< Chf^Wvo - QB(xo,r)Vo|||(xo,r-). 
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where we have applied the domain inverse inequality (see the Appendix in [21) 1 in the 
third line. Now applying Lemma [9.91 to the term ||vo — QB(xo,r)Vo|lB(xo r) obtain 


l|vo -QB(xo.r)Vo||iT < C'^T Ag ||£(vo) < C/lT ||e(vo) Hy. 

Combining (19.171) with the above estimate gives (19.161) . This completes the proof. □ 

Lemma 9.11. Assume that the finite element partition Th of is shape regular, 
and the finite element space Vh is constructed as in Section\^with fc > 1. Let 1 < r < 
k. Then, for any w G , p G and v G Vh, the following estimates 

hold true 

(9.18) |s((3/tw,v)| <C'/i’'||w||r+i|||v|||, 

(9.19) |C(v)| <C'/i’'||w||r+i|||v|||, 

(9.20) |0p(v)| <C'/inip||.|||v|||, 


where •^w(’) md Op(-) are given as in lOI) and 

Proof To derive (I9.18|) . we use (13.31) . the Cauchy-Schwarz inequality, and the 
estimates m and (19.31) to obtain 

|s((3hw,v)| =1 ^ hf^{QbQo'w - Qb'w,Qb'vo-Vb)aT 
T&Th 

= 1 ^ /iT^(QoW - W,(3bVo - V{,)aT 
T€Th 

-( f^T^WQb^O-^bfoT)^ 

T&Th TeTh 

<c( ^ /iT^IIQow - will + ||V(Qow - w)|||^ " |||v||| 

TgTh 

<Cft.’'||w||r+i|||v|||. 


To prove (19.191) . we use the Cauchy-Schwarz inequality, the estimates (|9.1D and 
(19.41) to get 


l^w(v)| 


“ Qhe(w))n)aT 

TdTh 


<c( XI - ^bWlr)^ { X ^'rlk(w) - Q,ie(w)|||j,) " 

TeTh TGTh 

||w||r+i( X ^T^lvo-Vblllr)"- 

Ten 


(9.21) 
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For the term ~ v&||gy, from Lemma 19.101 we have 

hT^\Wo-^b\\oT 

Ten 

<2 h;^^||vo - Q&voIIot + 2 Y - Vblla-r 

(9.22) Ten Ten 

<c Y II^K)||t + 2 Y ^T^WQb^O - ^bfaT 

Ten Ten 

<C|||vf. 

Substituting ()9.22p into ()9.21l) yields 

|^w(v)| < C'/l’'||w||r+l|||v|||. 


As to (19.201) . we use the Cauchy-Schwarz inequality, the estimates (|9.5I) and (I9.22p 
to obtain 


l^'pWI = 


Y (vo - Vb,(p- Qhp)n) 


OT 


Ten 


< 


( Y - Vblllr) " ( Y ^tUp-Q hpMlr)' 


Ten 

< Ch^\\p\\r\M\. 


Ten 


This completes the proof of the lemma. □ 


10. Numerical Results. In this section, we shall report some numerical re¬ 
sults for the weak Galerkin finite element scheme proposed and analyzed in previous 
sections. In our numerical investigation, we considered the linear elasticity equation 
dnj in the two dimensional square domain 17 = (0,1)^. The weak Galerkin approx¬ 
imations are obtained by using the lowest order finite elements. More precisely, the 
finite element space is given either by 

(10.1) 14 = {v = {vo,Vb} : Vo € [Pi{T)f,Vh G PRM{e), e C dT, T £%} 
or by 

(10.2) K = {v = {vo,v„} : Vo G [Pi(T)]2,v, G Pi(e), e C dT, T e %}■ 

Since PRM{e) C Pi(e) for each edge e C dT, then 14 is clearly a subspace of 14• It 
can be seen that all the theoretical results developed for 14 can be extended to 14 
without any difficulty. 

The discrete weak gradient and the discrete weak divergence are computed on 
each element T € Th according to their definition. More precisely, for any v G 14 or 
V G 14, the discrete weak gradient Vu,v on the element T is given by 

(10.3) (V„v, (/j)t = (vb,(/?- n)aT, V v? G [Po(P)]^''^- 
The discrete weak divergence V^, • v on T is given by 

(10.4) {Vyj ■v,iIj)t = {■Vb,'iljn)oT, 


V G Po(r). 
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Denote by U/j = {uo,U{,} and u the solution to the weak Galerkin formulation 
(1531) and the original equation (1131), respectively. Define the error by = QhU—Uh = 
{gq, et} where Q/jU is the projection of the exact solution u in 14 or 14, as 
appropriate. The error for the weak Galerkin finite element solution is computed in 
three norms defined as follows 


\\eh\ 


1= ^ / \Swieh)\'^dT + \ [ |Vu, • eh?dT + hj} f \Qbeo - efcpdsj , 

V. ^T" JT JdT J 


lleof = E / 


TeTh ■ 


l|eh|P = ^ hr j lebpds. 


TeTh 


dT 


It can be seen that |||e;i|||^ is an i/^-like norm for the error function, and ||eo|| and 
I jet, 11 are typical norms for the error in the interior and on the boundary of each 
element. 


10.1. Test Problem 1. Consider the elasticity equation dm) in the square 
domain 17 = (0,1)^ which is partitioned into uniform triangular mesh 7^ with mesh 
size h. The right-hand side function f is chosen so that the exact solution is given by 

_ f sin(a;) sin( 2 /) 

"= V 1 


Table 110.11 illustrates the computational result when the rigid motions are em¬ 
ployed on the boundary of each element; i.e., the space I 4 . Table [TU3] shows the 
result when linear functions are used on the boundary of each element; i.e., the space 
14- Theoretically, these two methods have the same order of convergence which is 
confirmed by these two tables. Note that the WG method with 14 has less number 
of degrees of freedom than that of t4. 

A numerical scheme is said to be convergent with order a if the error decreases 
proportionally to /i“, where h is the mesh parameter. From Tables lT0.il and [10.21 we 
see that the convergence of the weak Galerkin finite element scheme in the L^-norm 
is of order 2 and that in the iT^-norm is of order 1. The numerical results are in 
consistency with the theoretical prediction. 


Table 10.1 

Problem 1: X = I, /i = 0.5, and WG with 14 7^6 rigid motion space on element boundary). 


1/h 

||uo - Qou|| 

order 

||u{, - Qfeull 

order 

|||u/j - Qhu\\\^ 

order 

2 

0.0750 

- 

0.0424 

- 

0.3103 

- 

4 

0.0192 

1.97 

0.0115 

1.88 

0.1566 

0.99 

8 

0.0049 

1.98 

0.0031 

1.87 

0.0787 

0.99 

16 

0.0012 

1.99 

0.0008 

1.93 

0.0394 

1.00 

32 

0.0003 

2.00 

0.0002 

1.97 

0.0197 

1.00 




















Table 10.2 

Problem 1: X = I, n = 0.5, and WG with Vh (linear functions on element boundary). 
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1/h 

\\Uh - Qou|| 

order 

||uh - Qf,u|| 

order 

iUh - Q/tU||| 

order 

2 

0.0743 

- 

0.0424 

~ 

0.3082 

- 

4 

0.0190 

1.96 

0.0113 

1.90 

0.1555 

0.99 

8 

0.0048 

1.98 

0.0031 

1.88 

0.0782 

0.99 

16 

0.0012 

1.99 

0.0008 

1.93 

0.0392 

1.00 

32 

0.0003 

2.00 

0.0002 

1.97 

0.0196 

1.00 


10.2. Test Problem 2. In the second test, the linear elasticity equation (EU 
is also defined on the unit square domain 17 = (0,1)^, but the exact solution is given 
as follows 


u = 


(x + y)^ \ 
(x -yf ) ' 


The right-hand side function f is computed to match the exact solution. The numerical 
results, as shown in Tables ITO.31110.41 are again based on the uniform partition Th for 
this domain. The numerical results confirm the theoretical prediction developed in 
previous sections. 


Table 10.3 

Problem 2: X = 1, y, = 0.5, and WG with (the rigid motion space on element boundary). 


1/h 

||uo - Qoull 

order 

||u{, - Qbull 

order 

|||u/j - Q/tu|||^ 

order 

2 

0.4334 

-- 

0.2089 

- 

1.8410 

-- 

4 

0.1095 

1.98 

0.0528 

1.98 

0.9264 

0.99 

8 

0.0275 

2.00 

0.0133 

1.99 

0.4644 

1.00 

16 

0.0069 

2.00 

0.0034 

2.00 

0.2324 

1.00 

32 

0.0017 

2.00 

0.0008 

2.00 

0.1162 

1.00 


Table 10.4 

Problem 2: X = 1, p, = 0.5, and WG with Vh (linear functions on element boundary). 


1/h 

\\Uh - Qou|| 

order 

||uh - Qf,u|| 

order 

|||u?i - Q/jull 

order 

2 

0.4310 

- 

0.2009 

- 

1.8355 

~ 

4 

0.1088 

1.99 

0.0502 

2.00 

0.9233 

0.99 

8 

0.0273 

2.00 

0.0126 

2.00 

0.4628 

1.00 

16 

0.0068 

2.00 

0.0032 

1.99 

0.2316 

1.00 

32 

0.0017 

2.00 

0.0008 

2.00 

0.1158 

1.00 


10.3. Locking-Pree Tests. In the locking-free investigation, the linear elastic¬ 
ity equation (HH) has exact solutions given by 


/ sin(a;) sin(y) V ^-i / a; 
cos{x)cos{y) J \ y J 


The right-hand side function f is computed to match the exact solution (note that 
it is A-dependent). The numerical results are based on the same uniform partition 
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Th for the unit square domain. The results, as shown in Tables llO.SmO.T^ clearly 
indicate a locking-free convergence for the weak Galerkin finite element method in 
various norms, which is consistent with theory. 


Table 10.5 

WG with Vh based on the element of {Pi{T) / Pum {e)}, fJ, = 0.5, and A = 1. 


1/h 

\\uh - Qou|| 

order 

||ub - Qbu\\ 

order 

|||u/i - Q?iu||| 

order 

2 

0.0352 

- 

0.0331 

- 

0.1544 

- 

4 

0.0097 

1.86 

0.0120 

1.46 

0.0834 

0.89 

8 

0.0026 

1.91 

0.0037 

1.68 

0.0433 

0.94 

16 

0.0007 

1.96 

0.0010 

1.87 

0.0220 

0.98 

32 

0.0002 

1.98 

0.0003 

1.96 

0.0110 

0.99 


Table 10.6 

WG with Vh based on the element of {Pi{T)/PjiM{£)}, h = 0.5, and A = 100. 


1/h 

\\uh - Qou|| 

order 

||uh - Qfcull 

order 

|||u?i - Q/iu||| 

order 

2 

0.0344 

- 

0.0291 

- 

0.1449 

- 

4 

0.0100 

1.79 

0.0113 

1.36 

0.0774 

0.90 

8 

0.0028 

1.82 

0.0038 

1.59 

0.0403 

0.94 

16 

0.0008 

1.91 

0.0011 

1.81 

0.0205 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0103 

0.99 


Table 10.7 

WG with Vh based on the element of {Pi{T)/P[ im{s)}, h = 0.5, and A = 10,000. 


1/h 

\\uh - Qou|| 

order 

||ub - Qbu\\ 

order 

|||u?i - Q?iu||| 

order 

2 

0.0344 

- 

0.0290 

- 

0.1447 

- 

4 

0.0100 

1.79 

0.0113 

1.36 

0.0773 

0.90 

8 

0.0028 

1.82 

0.0038 

1.59 

0.0403 

0.94 

16 

0.0008 

1.90 

0.0011 

1.81 

0.0205 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0103 

0.99 


Table 10.8 

WG with Vh based on the element of {Pi{T)/PiiM{e)}, fi = 0.5, and A = 1,000,000. 


1/h 

\\Uh - Qou|| 

order 

||uh - Qbu\\ 

order 

|||u?i - Q/iu||| 

order 

2 

0.0344 

- 

0.0290 

- 

0.1447 

- 

4 

0.0100 

1.79 

0.0113 

1.36 

0.0773 

0.90 

8 

0.0028 

1.82 

0.0038 

1.59 

0.0403 

0.94 

16 

0.0008 

1.90 

0.0011 

1.81 

0.0205 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0103 

0.99 



































Table 10.9 

WG with Vh based on the element of {Pi(T)/Pi(e)}, n = 0.5, and A = 1. 
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l/h 

||u/i - Qou|| 

order 

||uh - Qhu|| 

order 

|||u/i - Q?iu||| 

order 

2 

0.0341 

- 

0.0313 

- 

0.1518 

- 

4 

0.0093 

1.87 

0.0115 

1.45 

0.0816 

0.90 

8 

0.0025 

1.91 

0.0036 

1.67 

0.0424 

0.95 

16 

0.0006 

1.96 

0.0010 

1.86 

0.0215 

0.98 

32 

0.0002 

1.98 

0.0003 

1.95 

0.0108 

0.99 


Table 10.10 

WG with Vh based on the element of {Pi(T)/Pi{e)}, = 0.5, and A = 100. 


l/h 

\\Uh - Qou|| 

order 

||uh - Qhu|| 

order 

|||u/i - Q/iu||| 

order 

2 

0.0340 

- 

0.0281 

- 

0.1441 

- 

4 

0.0098 

1.79 

0.0111 

1.34 

0.0769 

0.91 

8 

0.0028 

1.82 

0.0037 

1.58 

0.0400 

0.94 

16 

0.0007 

1.90 

0.0011 

1.81 

0.0204 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0102 

0.99 


Table 10.11 

WG with Vh based on the element of {Pi{T)/Pi{e)}, = 0.5, and A = 10, 000. 


l/h 

\\Uh - Qou|| 

order 

||uh - Qbu\\ 

order 

|||u?i - Q?iu||| 

order 

2 

0.0340 

- 

0.0280 

- 

0.1439 

- 

4 

0.0098 

1.79 

0.0111 

1.34 

0.0768 

0.91 

8 

0.0028 

1.82 

0.0037 

1.58 

0.0400 

0.94 

16 

0.0007 

1.90 

0.0011 

1.81 

0.0203 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0102 

0.99 


Table 10.12 

WG with Vh based on the element of {Pi (T)/Pi (e)}, ^ = 0.5, and A = 1,000,000. 


l/h 

\\Uh - Qou|l 

order 

||uh - Qhu|| 

order 

|||u;i - Q/iU||| 

order 

2 

0.0340 

- 

0.0280 

- 

0.1439 

- 

4 

0.0098 

1.79 

0.0111 

1.34 

0.0768 

0.91 

8 

0.0028 

1.82 

0.0037 

1.58 

0.0400 

0.94 

16 

0.0007 

1.90 

0.0011 

1.81 

0.0203 

0.97 

32 

0.0002 

1.96 

0.0003 

1.93 

0.0102 

0.99 
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